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1 Introduction 

The first purpose of this paper is to propose a formulation of general area-preserving 
motion of polygonal curves by using a system of ODEs. Solution polygonal curves 
CN ■ belong to a prescribed polygonal class, which is similar to admissible class used in the 

so-called crystalline curvature flow. Actually, if the initial curve is a convex polygon, 
then our polygonal flow is nothing but the crystalline curvature flow. However, if 
the initial polygon is not convex and does not belong to any admissible class, then 
the polygonal flow cannot be regarded as a crystalline curvature flow. Because the 
prescribed polygonal class is determined by the initial polygon and one can take any 
polygon as the initial data. On the other hand, in the framework of the crystalline 
curvature flow, the initial polygon should be taken from the admissible class. 

The second purpose is to discretize the ODEs implicitly in time keeping a 
given constant area speed, while the solution polygonal curve exists in the prescribed 
polygonal class. 

The organization of this paper is as follows. In the next section, we will 
introduce notion of polygonal motion and a polygonal class will be given. In section 
3, our problem will be formulated and some examples will be given. In the last 
■^j- ■ section, we will propose a scheme which achieves the second purpose and show 

convergence of the scheme. 

(N 
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C\l ■ 2 Polygons and polygonal motions 

o 

g§ ; 2.1 Polygons 

We define a set of polygons in R 2 : 

V := {T; T is a polygonal Jordan curve in M. 2 }. 

For r G V, the bounded interior polygonal domain surrounded by T is denoted 
by Q. For simplicity, we consider the case that Q is simply connected, but many 
of the following arguments are valid in other geometrical situations, even in three 
dimensional case, with some minor changes. 

Let r G V be an iV-polygon. The N vertices of T are denoted by Wj G M 2 
for j = 1,2, ... ,N counterclockwise, where w = and itfjv+i = Wi. Hereafter 
we use the periodic boundary condition F = Fn and Ftv+i = F x for any quantities 
defined on iV-polygon. 

The jth edge between Wj-i and Wj is 

T j = {(l-9)w j . l + 9w j ; 0<fl<l} (j = l,2,...,iV), 

and their lengths are |r^| := \iVj — We define Xj £ L°°(r) as 

'■""ft "r\r, <<-« ™ 



which is the characteristic function of Tj. 

The counterclockwise tangential unit vector and the outward unit normal 
are denoted by tj and rij, where tj = (wj — Wj^i)/\Tj\ and rij is defined such as 
det(n,-, tj) = 1. The outer angle at the vertex Wj is denoted by (fj G (— ir, it) \ {0}. 

We remark that 

cos (fj = tj + i ■ tj = rij + i ■ rij. 

We define the height of Tj from the origin hj := Wj ■ rij. They satisfy the equalities 

|r\j| = aj-xhj-i + bjhj + ajh j+1 (j — 1,2, ... , N), (2.1) 

where Oj := (sin Pj)^ 1 and bj := — cotipj-i — cot (fj. This equality can be checked 
from the fact that the straight line including r\, is expressed by the equation rij-x — 
hj. The total length of T is given by 

N N N 

l r l := Yl l^'l = J2( a i + b i + a i~i) h i = J^Vjhj, (2.2) 
j=i j=i j=i 

where rjj := cij + bj + Oj_i = tan(ipj/2) + tan(ipj_ 1 /2). 

The area of interior domain Q is denoted by which is given by 

1 N 
i=i 

The above symbols are also written as Q = fi(r), rij = rij(T) and hj = hj(T) 
etc., if we need to distinguish from quantities of the other polygons. 

2.2 Motion of polygons and Lipschitz mappings 

We consider a moving polygon T(t) e "P, where the parameter t (we call t time) 
belongs to an interval I C M. For 6 NU {0}, we call a moving polygon T(t) 
belongs to C fc -class on X, if the number of edges of T(t) does not change in time and 
Wj E C*(J, M 2 ) for all j = 1, 2, . . . , N. 

If fc > 1, we can define the normal velocity at x e Tj(t) which is the jth 
edge of T(t). We suppose x* e Tj(t*) and x* — (1 - 9)w j - 1 (t*) + 0wj(t*) for some 
G (0, 1), and define z:(#,t) := (1 - 9)w j - 1 (t) + 0Wj(t) G T(t). Then the outward 
normal velocity of Tj(t*) at x* is defined by 

Vj(x*,t*) := x(6,t*) ■ rij{t*) = (1 - fl)^.^**) • nj(t*) + ^(t*) • n^t*). 

Here and hereafter, the (partial) derivative of F with respect to t is denoted by F. 
We remark that Vj(-,t) is a linear function on each Tj(t). We define the normal 
velocity of F(t) by 

N 

n-,t) :=£^(,%(-,t)6r(r(t)), 

where Xj('it) ^ i s the characteristic function of Tj(t). 

For a C fc -class moving A^-polygon T(t) (t G X) with its interior domain Q(t), 
we construct Lipschitz mappings smoothly parametrized by t from a fixed polygonal 
domain. 



(2.3) 



We fix t* G X and define T* := T(t*) and fl* := We also choose another 

domain Q with \J te jQ(t) C Q an d fi x h. Our aim of this subsection is to construct 
Lipschitz mappings $(£) = $(•,£) from Q* to fi(£) smoothly parametrized by t. 



Proposition 2.1 Under the above condition, there exists e > 0, 1* := X D (f — 

£, i* + e), and $ G C k (l*, W X, °°(Q, M 2 )), and t/iey satisfy the following conditions. 

(1) zs a bi-Lipschitz transform from Q onto itself, i.e. $(£) zs bijective from 
Q onto itself and $(t) and are froia Lipschitz continuous on Q. 

(2) $(cc, t) = cc /or x in a neighborhood of dQ for t G X*. 



(3) <E>(f2*,£) = O(t) /or t G X* ; ana 1 $(£) is an affine map from each edge T* : = 
Tj(t*) onto Tj(t) with $(v)j(t*),t) = Wj(t). 

Proof Without loss of generality, we assume that Q is a bounded polygonal domain, 
too. We consider a triangulation T of Q* and Q. Namely T is a collection of 
triangular subatomics of Q with 

Q= \JK, and K n K' = if AT, AT' G T, K ^ K' , 

KeT 

and K K' is either the empty set , a common vertex or a common edge of K and 
K' G T, and there is a subset % C T such that 

= [JI, and r* n M = {w*}f =1 , 

A'eTo 

where Af denotes the set of all vertices of triangles in T and w* := Wj(t*). We also 
suppose that there does not exist any K G T with K n T* 7^ and .ff D <9Q 7^ 0. 
We assume that <&(x,t) has the following form: 

$(x,t) = A K (t) (fj (xeKET), (2.1) 

where Ax(t) is a 2 x 3 matrix depending on K G T and £. To determine A^(£), we 
suppose the condition: 

9(xt )-f x if^eAT\r*, 
na3,tj ~ if x = i»; GjVnr*. 

For sufficiently small e > 0, for t G X*, $(a;,£) is uniquely determined by the 
conditions (J21D and (Q. It is also clear that $(£) = $(•,£) G HA°°(Q ; R 2 ) is 
bijective from Q onto itself and $(x, t) = x for x G K if .K" fl T* = 0. 

Let us fix f G T (with AT fl T* 7^ 0) and let x\, x%, x 3 be the vertices of 
K. Then, from CT) . := $(scj,t) satisfies the condition y ; G C fe (X*,M 2 ) for 

/ = 1,2, 3. Since 



tt = ( (1 = 1,2,3), 



we have 



^(t) = (yi(t),?/ 2 (t),y3(t))( iC 1 1 ^ 



From this expression, we obtain A K £ C k (X*, M 2x3 ) and hence $ £ C k (T*, W 1 '°°(Q, 
follows. I 
As an application of this proposition, we have the following theorem which is 
well-known in the case T(t) is a smooth Jordan curve. 

Theorem 2.2 Let T(t) £ V be a C 1 -class moving N -polygon on an interval t £ X 
with its interior domain Q(t). For all £ C 1 (M 2 x 1), the map [t h- > (f>(x, t) dx] 
belongs to C 1 (X) and 



d_ 

dt 



I <p(x,t)dx= / <fi(x,t)dx+ / <f)(x,t)V(x,t) ds 
Jn(t) Jn(t) Jr(t) 



holds. 



Proof. We define f(t) := f a ^(j)(x,t) dx. Under the setting of Proposition 12.11 we 
show / £ C 1 {1*) and calculate /(**). For t el*, we have 

f(t)= [ <f>mx,t),t)J{x,t)dx, 
Jq* 

where J(x,t) is the Jacobian defined by J(x,t) := det(V x Q T (x, t)). We remark 
that J £ C\1\L°°(Q)) and j(x,t*) = diy x ^(x,t*) since $ £ C\l*, W l '°°(Q, R 2 )) 
and $(a;,t*) = a?. Hence, we obtain 

/(f) = I (V x <p(x,t*) ■ §{x,t*) +<j)(x,t*))dx+ I (f)(x,t*)j(x,f)dx 

div x ((f)(x,t*)®(x,t*))dx+ / <fr(x,t*)dx 

Jn* 



j)(x,t*)V(x,t*)ds+ / <f>(x,t*)dx 
i: in* 



In the case where = 1, we obtain the following formula for C 1 -class moving 
polygon particularly. 



d 'Q(t)\ = I V(x,t)ds, (2.3) 
r(t) 



where \£l(t)\ stands for the area of Q(t). 

We remark that Proposition 12.11 and Theorem 12.21 with their proofs are valid 
even in the three dimensional case. 



2.3 Polygonal motion 

For two polygons T and £ £ V, we define an equivalence relation T ~ S. We say 
r ~ S, if their numbers of edges are same (let it be N) and nj(T) = n>j(£) for 
all j = 1,2, ... ,N after choosing suitable counterclockwise numbering for T and E. 
The equivalence class of T £ V is denoted by V[T] := {£ £ V; £ ~ T}. 

We fix an A^-polygon T* e V and let V* := V[T*]. For T and S in V* , we 
define the distance between them by 

<i(r,S):= max \hi(T) — hAT,)\. 

j=l,2,...,N J 



Then, it is clear that (V*,d) becomes a metric space since it is isometrically em- 
bedded in WL N equipped with maximum norm | • |oo by the height function h: 
V* 3 T i-> h(T) = (h 1 (T),h 2 (T),...,h N (T)) T G R N . The following proposition 
is clear. 

Proposition 2.3 The set h(V*) is an open subset ofM. N . 

For any T and T 1 G V* and for 9 G [0, 1], we define 

h e := (1 - 9)h(T°) + ^(r 1 ) g R N . 

If there exists T e G V* with h(T 9 ) = h e , T e is called ^-interpolation of T° and T 1 . 
The ^-interpolation of T G V* and T 1 G V* is denoted by (l-^)r°+^r 1 := T e G V*. 
We remark that it satisfies 

|r?|>min{|r°|, |r}|}. 

For r G V* and e > 0, e-ball in = P[r] with center T is denoted by 
B(T,e) := {£ G V[T}; d(E,T) < e}. 
For an open set O C V* and r G 0, we define a positive number p(T, O) > as 

p(r,C) :=inf{|o-Mr)|oo; oeR JV \h(0)}. 

We remark that p(-, (9) is Lipschitz continuous with Lipschitz constant 1: 

\p(r,o)-p(z,o)\<d(r,z) (r, e g o). 

For a compact set JC C 0, we also define 
p(/C,0) :=minp(r,(9). 

Let a* := a,[r*] and 6* := 6j[r*]. Then, from the formula (12.11) . we obtain 

= lo^cvitr] - ViPl) + &KMr] - fyPD + s*(Vi[r] - VuPDI 
<c*rf(r,s) (j = i,2,...,iv), 

where we define 

CT:= i= ga v {|or_ 1 | + |6?| + |a?|}. (2.4) 

For a compact set JC C V* , we define 

<r(/C) :=min{|r.,|; T G /C, j = 1, 2, . . . , N} > 0. 

We consider a C fe -class moving polygon r(£) £ V* (t G X). We call it polyg- 
onal motion in "P* in this paper. We remark that a polygonal motion T(t) G "P* 
(t G X) belongs to C fc -class if and only if hj G C k {l) for j = 1, 2, . . . , N. If T(t) is 
a C 1 -class polygonal motion, its normal velocity Vj of T^t) is a constant on each 
Tj(t) and it is given by Vj(t) = hj(t). The formula (I2.3P is written in the form: 

d N 

Jt \m\ = J2\ T MvAt)- (2.5) 



We fix an equivalence class V* of polygons and let its j'th outward unit normal 
be rij and outer angle tpj. For r G V*, the polygonal curvature k,j of Tj is defined 
by 



We also define the polygonal curvature of T by 



N 
3=1 



The reason why this is called "curvature" is shown by the following proposition. 
Proposition 2.4 Let T(t) (t G X) be a C l -class polygonal motion in V* . Then 

d N r 

-\r(t)\ = J2\ T M^(m(t)= / <x,t)v(x,t)ds. 

dt j= i Jr® 



Proof. We obtain 

N N N 

ill' (It 



j=i j=i j=i 



from the formula (12.21) . 



3 Initial value problem of polygonal motion 
3.1 General polygonal motion problems 

We fix an equivalence class of iV-polygons V* as in §2.31 For an open set O C V* 
and T* G (0, oo], let F be a given continuous function from O x [0, T*) to R with 
the local Lipschitz property: For arbitrary compact set /C C and T G (0,T*), 
there exists L(/C,T) > such that 

|F(r, t) - F(E, t) U < L(K, T) d(T, S) (r, E G /C, t G [0, T]). (3.1) 

Under the condition (13.11) . for a compact set /C C and T G (0, T*), we also define 

M(/C,T) :=max{|F(r,t)| 0O ; T G /C, i G [0,T]} > 0. 

We consider the following initial value problem of polygonal motion. 

Problem 3.1 For a given N -polygon T* G O, find a C^^-class polygonal motion 
T(t) G O (0 < t < T < T») such that 

V j (t) = F j (T(t),t) (te[0,T\, j = l,2,...,N) 

r(o) = r*. 



Under the Lipschitz condition (I3.ip . it is clear that there exists a local solution T(t) 
in a short time interval [0, T], since Problem 13. II can be expressed by an initial value 
problem of an ordinary differential equations for h(t) = (hi(t), h,2(t), . . . , hiy(t)) T . 
We often assume the following condition for Ff. 

N 

^ir^r,*)^ (TeO, te[o,T„)), (3.2) 

i=i 

where \i is a fixed real number. Under the assumption (13.21) . from the formula (12.51) . 
any solution T(t) to Problem 13.11 has the property of the following constant area 
speed (CAS in short): 

The polygonal flow is regarded as the crystalline curvature flow if the initial 
polygon T* is convex (as mentioned in introduction, if T* is not convex, then the 
polygonal flow is different from the crystalline curvature flow). There are many 
articles about the crystalline curvature flow and asymptotic behavior of solutions 
[H El [H [51 El [91 HOI [Til [T2], etc., which started from the pioneer works [2] and [7] . 

3.2 Examples of problems of polygonal motion 

In this section, we give some basic examples of polygonal motions which are nice 
polygonal analogues of corresponding smooth moving boundary problems. 

Problem 3.2 (polygonal curvature flow) For a given N -polygon T* G V* , find 
a C 1 - class family of N -polygons U 0<t<T F(t) c (T < T*) satisfying 

V j {t) = -K j (t) (te[0,T], j = 1,2,..., TV), 

r(o) = r*. 

The solution has CAS property with a = —2^2^ =l t&n((pj/2): 

d N N N 

j t \m\ = -5>j(*)|Ii(t)| = = -2^tan^ = const. 

3=1 i =1 i =1 

Problem 3.3 (area-preserving polygonal curvature flow) For a given N- 
polygon T* G V* , find a C 1 -class family of N -polygons Uo<t<T-' n (^) C ?* (T < T*) 
satisfying 

'Vj(t) = («(-, t)) - Kj(t) (t G [0, T], j = 1, 2, . . . , AT), 

r(o) = r*. 

Here («(•, t)) is the mean value of k on T(i): 

( K (.t))-—[ rJx t)dc- ^= lVi - 2 E t =itan(^/2) 



The solution has CAS property with \x = 0: 

±Mt)\ = (K(;t))\T(t)\- [ K( X ,t)ds = 0. 

at Jrtt) 



r(t) 

In what follows, the mean value of F on the edge Tj is denoted by 



(F),:=-L ! f( x )ds. 

I .7 1 J r\ 



j\ JY j 

Let G be a bounded Lipschitz domain in M 2 . We define 
O g := {T G V*; n(T) D G}. 

Problem 3.4 (area-preserving polygonal advected flow) Let us consider a di- 
vergence free vector field u G C^M 2 \ G; R 2 ) with divu = in R 2 \ G. For a 
given N -polygon T* G Oq, find a C l -class family of N -polygons Uo<t<T-^(^) c ®G 
(T < T*) satisfying 

'V j (t) = {u-n) j (te[0,T\, j = l,2,...,N), 

r(o) = r*. 



The solution has CAS property with \x — n ■ uds: 

IdG 



d 



\fl(t) | = — / rij ■ uds = — div udx + n ■ uds = n ■ uds, 

J r(t) Jn(t)\G JdG JdG 



dt Jv(t) Jn(t)\c 

where rij and n are the outward unit normal vector to d(Q(t) \ G). 

Problem 3.5 (polygonal Hele-Shaw flow) For a given N -polygon T* G Og and 
a function b defined on dGx[0, T], find a C 1 -class family of N -polygons IJo<t<T ^(0 C 
O g (T < T*) satisfying 

, V j (t) = -n j -Vp(x,t) (xeFjit), te[0,T}, 3 = 1,2,. ..,N), 
A P (x, t) = (xe Q(t) \G, te [0, T]), 

(p(-,*)>j = K iW (te[0,T}, j = 1,2,..., AO, 

J^(a;,f) = &(*,*) (zGdG, te[0,T]), 

r(o) = r*. 

Here dp/dn = V x p-n and rij and n are the outward unit normal vector to d(Q(t) \ 
G). The solution has a given area speed property: 

^-\n(t)\ = - [ n r Vpds = - [ _Apdx+ [ ^-(x,t)ds= [ b(x,t)ds. 

dt Jr(t) JSl(t)\G JdG V n JdG 



If b(x, t) = bo for a given constant bo, then the solution has CAS property with 
pi = \dG\bo- 



4 Numerical schemes 



4.1 Notation 

In §H we consider time discretization of Problem 13.11 with the following notation. 
The discrete time steps are denoted by = t < t\ < t% < ■ ■ ■ < t r - n < T. The step 
size which may be nonuniform and their maximum size are defined by 

r m ■= Wi -t m (m = 0, 1, • - • , m - 1), r := max r m . 

0<m<m 

Approximate solution of T(t m ) is denoted by T m G V*. Quantities of the polygon V 
are denoted by hj := hj{T m ), and := Kj{T m ), etc. We define ef := hj(t m ) - h 
and e m := {e?, e%, . . . , e^) T G R*. Then we have d(T(t m ), T m ) = |e m |oo. 

The discrete normal velocity V™, which is an approximation of Vj(t 
hj(t m ), is defined by 

hf +1 - K 

Corresponding to the formula (12.51) . the following formula holds 



in 



ym. = _j_ L_ ( m = 0,1,..., m-1) 



|nm+l| lorn J* \V m \ -X- |r m+1 l 

1 - = E — 2 y r- (") 

• m ■ i ^ 

3=1 

This has a form of sum of areas of iV trapezoids and is actually derived from ( 12.3ft 
as follows: 

1 - 

3=1 



1 " 

- ^ {(|r; n+1 | + |r™ \){h™ +1 - hf) + \v™ +1 \hj - \vj\h 



2 

3=1 



m I tm+1 ' 
3 l Ai ? 



AT „ JV 



= y E (|r T +l1 + l r 3 m l) y 3 m + \ll (\ T T +1 K - \ r TK +1 ) ■ 

3=1 3=1 

where the last sum is equal to zero due to the equality (12.11) . 

In the following subsections, we suppose that there exists a unique solution 
r(t) for < t < T < to Problem 13.11 under the condition (13. ip , and that discrete 
time steps = t < t\ < t 2 < ■ ■ ■ < t m < T are given a priori such as the uniform 
time stepping t m = mr. It is, however, possible to apply any a posteriori adaptive 
time step control scheme. 

4.2 Euler scheme 

We consider the following Euler scheme to discretized Problem 13.11 

Problem 4.1 For a given N -polygon G O and time steps = t < tx < t 2 < 
' ' ' < tm < T , find polygons T m G O (m — 1,2, ... , fh) such that 

VJ n = F j (T m ,t m ) (m = 0,1,2,..., m-1, j = 1, 2, . . . , N), 

r° = r*. 



Theorem 4.2 We suppose the condition 113.1]) and that {r(t)}o<t<r be a C k+1 -class 
solution of Problem 13.11 for k = or 1. There exists 5* > 0, r* > 0, C > and a 
non-decreasing function u(a) > with 

[oil) ifk = 0, 

U ^ = \nf\ ti 1 0Sai °' (4 ' 2) 
[0(a) %fk = \, 

such that, ifd(T*, T°) < 5* and r < t* , then T m G O (m = 1,2,..., fn) is determined 
by the Euler scheme ( Problem ]^. 1\) and satisfies the estimate 

max d(T(t m ), T m ) < u{r) + Cd(T(0), T°). 

0<m<m 

The proof is similar to the one of Theorem 14.71 

4.3 Second order implicit scheme 

We consider the following implicit scheme for Problem 13. 1[ 

Problem 4.3 For a given N -polygon r, 6 and time steps = to < ti < h < 
■ ■ ■ < tfn < T, find polygons T m G O (m = 1,2,..., fn) such that 

V™ = F 3 (T m+1 / 2 , W1/2) (m = 0, 1, 2, . . . , fn - 1, j = 1, 2, . . . , N), 

r° = r„ 

where r m+1//2 and t m +i/2 are the 1 / 2-interpolations: 

r m+1/2 :=^-^ — ev*, Wi/ 2 := m 2 m+1 = *m + y- 

This is a generalized version of [8] for area-preserving crystalline curvature flow. 

Proposition 4.4 We suppose the constant speed condition f)3.2p . Let T m G O (m = 
1,2, ... ,fh) be a solution of Problem \A.?>\ Then it satisfies 

\Q m+1 \ = \Q m \+ fiT m (m = 0,1,..., m- 1). 

In other words, \Q m \ = \Q(t m )\ holds if the exact solution Q(t) of Problem \'S.l\ exists. 

Proof. Since |r™ +1/2 | = {\Tf\ + |r™ +1 |)/2, we have 

|O m+1 | - \tt m \ v 



E, r m+l/2|j ? / rm+ l/2 , \ _ 

3=1 

from the formula fl4.ll) . I 
Since Problem 14.31 is an implicit scheme, it is not clear whether r m+1 G O 
can be determined uniquely from the previous polygon T m G O, the time t m , and 
the time step size r m . Another question is how to solve the equations 

(pm 1 rtn+1 \ 
- 2 ,Wi/2 , (4-3) 

to obtain (approximation of) r m+1 numerically. 



We fix f G O and t G [0, T), which correspond to T m and t m +i/2, respectively. 
Let /C be a compact convex set in V* with f 6 K C 0. For S G /C and f G 
((^(f^M^T)" 1 ), we can define S G £> by 



h(t) = h(t)+tF 



s + r 



t . 



In other words, we can define A(X3) := A(E; f , t, f) := S which is a map from /C to 
(9. We have the following lemma. 

Lemma 4.5 Let e G (0, p(T, O)) and A G (0,1) be fixed. Suppose that f satisfies 
the condition 



< t < min <T-t, 



2A 



M(/C,T) L(/C,T) 

where JC := £?(r, e). T/ien A maps /C mto /C and 
rf(A(S 1 ), A(£ 2 )) < Ad(£\ S 2 ) (S 1 , S 2 G /C). 



(4.4) 



Namely, A is a contraction mapping on fC and there exists a unique fixed point of A 
in fC. 

Proof. It is enough to show (|4.4p . which is proved as follows: 
^(E 1 )^ 2 )) 



|MA(s 1 ))-MA(s 2 ))U = f 



s 1 + r 



s 2 + r 



< TL(K,T)d ( 
= fL(t,T) 



(y 1 



s x + r E 2 + r 

2 ' 2~~ 
/^I! 1 ) + /i(f ) h(£ 2 ) + h(f) 



= ^L(JC, T)d (E\ S 2 ) < Xd (S 1 , S 2 ) . 

I 

We immediately have the following theorem, which gives us efficient numerical 
scheme to obtain r m+1 . 

Theorem 4.6 Let K be a compact set in O and let e G (0, p(/C, O)). We define 



Sg/C 

For fixed m (< fh) in Problem \A.3\ we assume that T m G /C and 

e 2\ 

r m < mm 



M(/Q,T)'L(/Q,T) 



where A G (0, 1). Then there exists uniquely T m+1 G B(T m ,e) satisfying (14. 3p . 

Furthermore, r m+1 is i/ie /ised pomt o/ t/ie contraction A m := A(- ; T m , t m , r m ) 
in B(T m ,e), and is given by the limit of A m (T m ) as v — > oo wife i/ie following 
estimate 



d(T m+ \ A u m (T m )) < X u d(T m+ \ T m ) (i/GN). 



Theorem 4.7 We suppose that {r(t)} <t<T be a C k+l -class solution of Problem XS.W 
fork = 0, 1, or 2. There exists 5* > 0, t* > 0, C > and a non- decreasing function 
(jj(a) > with 



uj(a) 



o(a k ) if k = or I, 
0(a 2 ) ifk = 2, 



as a 10, 



(4.5) 



such that, if d(T*,T°) < 5* and r < r* , then T rn G O (m = 1,2,..., rh) are 
inductively determined and 

max d(T(t m ), T m ) < u(t) + Cd(T(0), T°), 

0<m<in 

holds. 

Proof. We put p := p({T{t); 0<t< T},0), and fix 5 G (0,p) and e G (0,p - 5). 
We define 



K:= |J fl(T(f),<5), /C e : = 



0<t<T 



se/c 



L:=L(JC E ,T), R(a) := e* fl - ^ 



-l/o 



(0 < a < 2/L), 



Pr, 



(m = 0, 1, 2, . . . , m), 



where a non- decreasing function a; (a) (0 < a < T), which satisfies ( 14.51) . will be 
defined in (14.91) later. Since R(-) is an increasing function and lim a j_ R(a) = e L , 
there exists 5* > and r* > such that 



R(t*) 5* + 



< S, t* < min 



M{K e ,Ty L 



For m = 0,1,2, ... ,m — 1, we will prove the following inductive conditions: 

T m G/C, Pm <R{r) tm P o =4- 3 r m+1 G/C, p m+1 < R(r)^po. (4.6) 

The conditions T° G /C and po < i?(r)°po for the case m = are obviously satisfied. 

Let us assume the conditions T m G /C and p m < R(r) trn p for a fixed m. 
Then, from Theorem 14.61 there exists r m+1 uniquely in B(T m ,e) C K, £ , and we have 

e m+1 - e m = h(t m+1 ) - h(t m ) - r m V m = r m {C + (V(t m+1/2 ) - V m ) } , (4.7) 
h(t m+ x) — h(t m ) 



- h{t 



m+l/2), 



where V(t) := (V^t), V 2 (t), V N (t)) T and V m := (Vf, V 2 m ,...,V%j . 

The last term of (14. 7p is estimated as follows. Since r(i m+ i/ 2 ) G /C C JC £ and 
T m+i/2 e _B(r m ,e) C K e , we have 



|V(t. 



m+l/2, 



- V* 



|F(r(t 



+1/2 J, t m +X/2 



-F(T m+l ' 2 ,t m+1/2 ) 



<Ld{T{t m+1/2 ),T m+1 ' 2 ) 



h(t 



h m + h m+1 



m+l/2) 



-(e m + e m+1 ) 



where 

. m h(t m ) + h(t m+1 ) 

Q '■= 2 n(t m+1 / 2 ). 

Combining (14.71) and (14.81) . we obtain 

' 1 



e m+1 l < le m l +r l£ m l +r L 

° |oo _ |° |oo 'm |S loo 1 m J -' 



Tm.L 



-(e m + e m+1 ) - C m 



^ I . t^t. I I m /|^™+1| i I _m | \ i ^ Html i T \ /-mx 

< l e loo + (|e |oo + |e |ooJ+^m(|€ |oo + -^|C I 



looy 



By the Taylor expansion, we can obtain an non- decreasing function cu(a) (0 < a < T) 
which satisfies the condition (14.51) and the inequality 

|rioo + £|C m |oo<^(r). (4.9) 

Hence, we have 

M |e m+1 L < (\ + Ic-lco + r r Mr) 



and this inequality is equivalent to 



Tm.L \ ( TffiL 



Pm+1 < 1 + -—- V 



From the inequalities 



T m lj \ I _ TLj \ , / _ T m lj \ Tm L/2 



'-Irjn'-Tj md ( 1 + ^J^ 

we obtain 

p m+ i <U- ^) 1 (l + ^) p m < i2(r) r »(i2(r)*-'pb) < R(r)^ Po . 
The condition r m+1 e /C follows from this estimate as 

|e m+1 U < Pm+1 < R(t*Y^ po < R(t*) T (V + ^P-) < 5. 



L 

Hence, we have proved (14.61) . which leads us to the estimate: 

\e |oo < H[t ) Me loo H — — I — < H(t ) |e ]oo H ^(r). 

The assertion of the theorem is obtained by putting C := R(t*) t and denoting the 
above term L~ 1 (R(r*) T — 1)u;(t) again by oj{t). I 
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